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Abstract 



A general analysis of squeezing transformations for two mode systems is given 
based on the four dimensional real symplectic group 5|?(4, K). Within the 
framework of the unitary metaplectic representation of this group, a distinc- 
tion between compact photon number conserving and noncompact photon 
number nonconserving squeezing transformations is made. We exploit the 
Sp{4:,^) — 50(3,2) local isomorphism and the U{2) invariant squeezing cri- 
terion to divide the set of all squeezing transformations into a two parameter 
family of distinct equivalence classes with representative elements chosen for 
each class. Familiar two mode squeezing transformations in the literature are 
recognized in our framework and seen to form a set of measure zero. Examples 
of squeezed coherent and thermal states are worked out. The need to extend 
the heterodyne detection scheme to encompass all of U (2) is emphasized, and 
known experimental situations where all U (2) elements can be reproduced are 
briefly described. 
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I. INTRODUCTION 



The theoretical analysis |T| and experimental 0, 0], realization of squeezed states of 
radiation continue to receive a great deal of attention. While much of the work so far has 
concerned itself with single mode situations [Q, p some analysis of two-mode states has also 
been presented 0, @]. Other nonclassical effects of radiation beyond second order have also 
received attention in the literature 0. More recently, a general invariant squeezing criterion 



for 77,-mode systems has been developed by some of us elsewhere [ITU 



The purpose of the present paper is to study squeezing transformations for two-mode 
systems, and to develop a classification scheme for them motivated by the above men- 
tioned invariant squeezing criterion. Basic to all such discussions is the four-dimensional 
real symplectic group Sp{4,^), of real linear homogeneous canonical transformations, and 
the unitary metaplectic representation of this group acting on the Hilbert space of states 
for a two-mode quantum system. The structure of the noncompact group Sp{4:, 3?) ( and 
its n-mode counterpart Sp{2n, 3?)) leads to a natural separation of its elements into passive 
(compact), and active (noncompact) types. Here the adjectives passive and active mean 
total photon number conserving and nonconserving respectively. This group theoretical 
framework gives us an unambiguous way of defining precisely the family of squeezing trans- 
formations; they are the active elements of 5*^(4, 3?) and they do not form a subgroup. The 
action of the maximal compact (passive) subgroup U (2) of 5*^(4, 3?) on the set of squeezing 
transformations by conjugation leads to a natural equivalence relation, leading to the emer- 
gence of equivalence classes and convenient representative elements as well. In studying the 
physical properties of a state subjected to squeezing, therefore, we are able to isolate the 
dependence on intrinsic squeezing parameters and separate them from other passive factors. 
As might be expected, the single squeeze factor encountered in the studies of single mode 
states gets enlarged here to two independent intrinsic squeeze factors; and it turns out that 
the two mode squeezing transformations so far studied in the literature form a very small 
subset of all the independent available possibilities. 
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The material in this paper is arranged as follows: Section |T| sets up the basic kine- 
matics for two mode systems, the Fock and Schrodinger representations and the actions of 
Sp{A, 5R) on the canonical variables and states. The variance matrix for a general state and 
its change under 5*^(4, 3?) are derived. After identifying the maximal compact or passive 
U{2) subgroup of 5^(4,3?), the f/(2) invariant squeezing criterion for two mode systems 
is discussed. Section |TTl| introduces the generators for the metaplectic representation of 
5*^(4,3?) and brings out the connection to the SO{3,2) Lie algebra. The photon number 
conserving compact generators and the remaining noncompact ones are clearly identified. 
Based on the polar decomposition theorem for general elements of 5*^(4, 3?) we are then led 
to a precise definition of squeezing transformations: These are single exponentials of linear 
combinations of the noncompact metaplectic generators. We then proceed to break up the 
set of all squeezing transformations into equivalence classes under U{2) action. We find that 
these classes form a continuous two parameter family describable by points in an octant in 
a two dimensional plane. The families of Caves-Schumaker transformations and essentially 
single mode transformations correspond to one dimensional lines bounding this octant and 
so are of measure zero. Section |rV| applies our formalism to two mode squeezed coherent 
and thermal states. In Section [V| we see how the heterodyne detection scheme fits into our 
analysis. We argue that it is necessary to experimentally realise all elements of the U{2) 
subgroup of 5*^(4, 3ft); the heterodyne scheme only handles a one parameter subset of U{2). 
Two examples of two-mode situations where all elements of U{2) can be experimentally 



realised are briefly described. Section VI contains some concluding remarks. 



II. SYMPLECTIC GROUP FOR TWO MODES AND THE SQUEEZING 

CRITERION 

We consider two orthogonal modes of the radiation field, with annihilation operators 
aj,j = 1,2, and corresponding creation operators Oj. These two modes could, for example, 
be two different frequencies for the same or different propagation directions and polarizations, 
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two different propagation directions at a common frequency, two different polarization states 
of plane waves degenerate in frequency and direction of propagation, etc. We arrange these 
operators in the form of a four-component column vector: 



r(c) 



(a' 



(c)^ 



ai 

02 



1,2,3,4. 



(2.1) 



The superscript (c) indicates that the entries here are complex, i.e., nonhermitian, operators. 
For discussion of quadrature squeezing, however, we need to also deal with the hermitian 
quadrature components of these operators. Therefore we define another column vector ^ 
with four hermitian entries, related to ^^'^^ by a fixed numerical matrix, as follows: 



Pi 

yP2j 



1 —i 



1 

72 



1 i 

1 i 

1 -i 
1 -i 



(2.2) 



The canonical commutation relations can now be written either in terms of ^ or (in two 
ways) in terms of ^('^^: 



/ 



(A 



'ab) 



10 
1 
-10 
0-100 



(2.3) 



'^ 1 ^ 
10 
0-10 
-1 

A general real linear homogeneous transformation on the g's and p's is described by a 
4x4 real matrix S acting as follows: 



/ 



(2.4) 



If the are to satisfy the same commutation relations as the ^a, the condition on S is: 



S (3S^ = f3. 



(2.5) 



This is the defining property for the elements of the group S'p(4, 3?), which is a noncompact 
group: 



Sp{4, ; 



{5 = 4x4 real matrix \ S (3 S'^ ^ (3} . 



(2.6) 



Note that P — — /3 ^ itself is an element of 5'p(4, 3?) whereas E3 is not. Further, S e -5^(4, 3?) 
imphes -S, S^, S-^ = (3 S^(3'^ e Sp{4, 3?) and det 5 = 1 for every S e Sp{4, 3?). 

The action of 5*^(4, 3fJ) on the nonhermitian operators is by a (generally) complex 
matrix 5'*^'^) related to 5" by conjugation with Jl: 

SeSp{A,U): ^'^S^^ 

^(c)/ ^ ^(c) ^(c)^ 

^nsnl (2.7) 
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We denote the Hilbert space on which ^ and ^'^'^^ act by Ti. Since the hermiticity prop- 
erties and commutation relations of the E,a are maintained by the transformation (|2.4| ) for 
any S G 5*^(4, 3?), and since the S,a act irreducibly on 7i, it follows from the Stone- von 
Neumann theorem that it should be possible to construct a unitary operator U{S) on 



Ti implementing ( p.4|) via conjugation: 

S E Spii, 3?) : SabCb = U{S)-^ U{S) , 

U{Sy U{S) = 10X171. (2.8) 

This U{S), for each S, is clearly arbitrary upto (and only upto by virtue of irreducibility) 
an S'-dependent phase. Even after making use of this freedom, however, it turns out that 
we cannot choose the operators U{S) for various S so as to give us a true unitary represen- 
tation of 5*^(4, 3ft) on Ti. Rather, after maximum simplification, they give us a two-valued 



representation of this group [12 



Si,S2e SpiA, 3ft) : UiSi) U{S2) = ±U{Si S2). (2.9) 

Alternatively we can regard the U{S), chosen so as to obey this composition law, as providing 
a true and faithful unitary representation of the four-dimensional metaplectic group Mp{4), 
which in turn is a two-fold cover of Sp{4, 3ft). In the literature the operators V({S) are often 
said to provide the metaplectic representation of 5*^(4, 3ft). The connection between Mp{4) 
and 5^(4, 3ft) has some similarity with the one between SU (2) and S0{3) familiar in angular 
momentum theory. 

Now we consider physical states of the two-mode system, the action of 5*^(4, 3ft) on them, 
and the statement of a suitable squeezing criterion. Let p be the density operator of any 
(pure or mixed) state of the two-mode radiation field. With no loss of generality we may 
assume that the means Tr (p^a) of vanish in this state. (Any non-zero values for these 
means can always be reinstated by a suitable phase space displacement which has no effect 
on the squeezing properties). Squeezing involves the set of all second order noise moments 
of the quadrature operators qj and pj. To handle them collectively we define the variance 
or noise matrix V for the state p as follows: 



V={V, 



ab 



(2.10) 



This definition is valid for a system with any number of modes. For a two-mode system it 
can be written exphcitly in terms of qj and pj as: 



/ 



V 



\{{(ii^Pi]) {(I2P1) 



{(I2P1) i({g2,P2}) 
{pi) {pm) 



\ 



\ 



(2.11) 



{(I1P2) |({g2,P2}) {P1P2) {pD 

This matrix is real symmetric positive definite and obeys additional inequalities expressing 



the Heisenberg uncertainty principle of quantum mechanics [IC]. 

When the state p is transformed to a new state p' by the unitary operator U{S) for some 
S G 5*^(4, 9ft), we see easily from eqs ( P78| , |2.1CI| ) that the variance matrix V undergoes a 
symmetric symplectic transformation: 



S e Sp{A, 3ft) : p' = U{S) p U{S)-^ 
V'= SV . 



(2.12) 



This transformation law for V preserves all the properties mentioned at eq. ( p.ll| ). 

Towards setting up a squeezing criterion, we identify an important subgroup of S'|j(4, 3ft), 
namely its maximal compact subgroup K = U{2). This consists of matrices S G 5*^(4, 3ft) 
having a specific block form determined by two-dimensional unitary matrices belonging to 
f/(2): 



K = U{2) = [S{X, Y) G Sp{4, ^)\U = X - iY e U{2)} C Sp{4, K), 



S{X,Y) 



^ X y\ 



-Y X 



(2.13) 



]Here X and Y are the real and imaginary parts of U; and the unitary condition W U = 1 
on U ensures that S{X,Y) obeys the symplectic condition ( p^) . We also recognise that 



this subgroup K is the intersection of the symplectic and the orthogonal groups in four real 
dimensions, namely 

K = U{2) = Sp{A,^)nO{A) : 

S{X,Y) {(3 oil) S{X,Yf = (3 or 1. (2.14) 

The complex form S^^\X, Y) corresponding to S{X, Y), and which acts on is rather 
simple; it is given by 

S^'\X,Y) = S^'^U) 

( X-iY ^ 

(2.15) 



y X + iY J 

Thus the maximal compact subgroup of S'p(4, 3fJ) mixes ai and 02 unitarily, but does not mix 
the Oj with the Oj. This fact may be explained by saying that the subgroup U{2) of Sp{4:, 3?) 
consists of passive, or total photon number conserving, transformations. 

We see that the physical requirement that the total number of photons be conserved 
singles out a unique maximal compact subgroup of 5*^(4, R) out of many equivalent ones. In 
contrast, elements of S'p(4, 3?) outside the subgroup U (2) are non-compact elements which do 
not conserve total photon number, and so describe active transformations. (These properties 
of compact and non-compact elements of 5*^(4, 3?) will become transparent when we identify 
their generators in the following Section). 



As has been discussed in detail elsewhere |T0[, H^l) ^ multimode system it is physically 



reasonable to set up a definition of squeezing which is invariant under the subgroup of 
passive transformations of the full symplectic group. For the present case of two- mode 
systems, we evidently need a f/(2)-invariant squeezing criterion. That is, our definition 
must be such that if a state p with variance matrix V is found to be squeezed, then the state 
U{S{X,Y)) p U{S{X,Y)y^ with variance matrix V = S{X,Y) V S{X,Yf must also be 
squeezed, for any U = X — iY E U{2). 

Conventionally a state is said to be squeezed if any one of the diagonal elements of V is 
less than 1/2. The diagonal elements correspond, of course, to fluctuations in the "chosen" 



set of quadrature components of the system. The [/(2)-invariant definition is as follows: 
the state p is a quadrature squeezed state if either some diagonal element of V is less than 
1/2 (and then we say that the state is manifestly squeezed), or some diagonal element of 
V = S{X, Y) V S{X, Y)^ for some U = X - lY e U{2) is less than 1/2: 

p is a squeezed state -v^ 
{S(X,Y)VS(X,Yf)^^<\ 

for some a and some X — iY G U{2). (2.16) 

That is, running over S{X, Y) & K is the same as running over all possible sets of quadrature 
components. We may say that since any element of U{2) passively mixes the two modes, 
the appropriate S{X, Y) & K which achieves the above inequality for some a (assuming the 
given V permits the same) just chooses the right combination of quadratures to make the 
otherwise possibly hidden squeezing manifest. 

To implement this definition in practice, it would appear that even if a state is intrinsi- 
cally squeezed, we may have to explicitly find a suitable U{2) transformation which when 
applied to V makes the squeezing manifest. This however could be complicated. Here 
the point to be noticed and appreciated is that diagonalisation of a noise matrix V gen- 
erally requires a real orthogonal transformation belonging to SO (4) which may not lie in 
U (2) = 0(4) n 5*^(4, 3fJ). It is therefore remarkable that, as shown in [|l^, the f/ (2) -invariant 
squeezing criterion ( p.l6| ) can be expressed in terms of the spectrum of eigenvalues of V, 
namely: 

p is a squeezed state 4^ i(y) = least eigenvalue ofV<-. (2-17) 

That is, while the diagonalization of V is in general not possible within K = U{2) which 
is a proper subgroup of 0(4), any one particular (and hence the smallest) eigenvalue of V 
can be made to become one of the diagonal elements of V transformed by an appropri- 
ate S{X, Y) G K. In other words any quadrature component can be taken to any other 
quadrature component by a suitable element of U{2). We shall in the sequel work with the 

10 



f/(2)-mvariant squeezing criterion (|2.16 



2T7D . 



III. THE 5P(4,5R) - SO{3,2) CONNECTION AND CLASSIFICATION OF 
TWO-MODE SQUEEZING TRANSFORMATIONS 

We have shown in the previous Section that the group Sp{4, 3?) of hnear canonical 
transformations contains two kinds of elements: passive total photon-number conserving 
elements belonging to the maximal compact subgroup K = U{2); and active noncompact 
elements lying outside this subgroup, and which do not conserve total photon number. It is 
clear from the f/(2)-invariant squeezing criterion ( p.l6| , |2rT7D that the former elements can- 



not produce squeezing. This is because the corresponding changes in the variance matrix, 
V = S{X, Y) V S{X, Y)'^ , being similarity transformations preserve the eigenvalue spec- 
trum of V] hence i{V) > 1/2 implies i{V') > 1/2 and conversely for every S{X,Y) G U{2). 
The noncompact elements of 5*^(4, 3?), on the other hand, while they do not form a subgroup, 
have the potential to produce a squeezed state starting from a nonsqueezed state. Thus they 
may be called squeezing transformations. The following questions then naturally arise: what 
are the really distinct squeezing transformations which are not related to each other by just 
passive transformations, and how can they be invariantly labelled or parametrised? 

To answer these questions it is useful to work at the level of the Lie algebra 5*^(4, 3?) of 
Sp{4, 3fJ), and the hermitian generators of the unitary metaplectic representation U{S) acting 
on 7i. As is well known, these generators are basically all possible hermitian symmetric 
quadratic expressions in the canonical variables q and p. They may be expressed more 
transparently for our purposes in terms of a's and a^'s. While it is possible to set up the 
generators in a uniform manner for all Sp{2n,^), in the two-mode case we can exploit the 
fact that the group 5^(4,3?) is locally isomorphic to the de Sitter group 50(3,2) [|1^], as 
they share the same Lie algebra. We therefore choose the basis for Sp{A, 3?) to make this 



explicit; it helps us visualise geometrically the analysis to follow |T5[. We define the ten 
generators Q, Jr, Kr, L^, r = 1, 2, 3 of S'p(4, R) as: 
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D — 

V — 


^ ( AT -L ^ \ — ^ ( n'^ n -L n'^ n -LiV 


l^O. id j 


7 


-\a1a2 + a2ai), 




7 


^ ^ t t ^ 




T 

-'3 — 


2 v^i^i 0202^, 


{6.1^} 


Ki — 


-(Oi +01-02 - 02j, 




A2 — 


V/7t2 „2 I _t2 2n 
-01 + 02 -02j, 




K — 


1 / f f N 

2 '^12 + Oi02j, 




7" 


*/'^t2 „2 „t2, „2\ 

-(Oi -Oi -02 +02), 






-'-/ t2 1 2 1 t2 1 2\ 
-(ol +0i + 0^ +O2), 




i^3 = 


^ ^ t t ^ 
-2('^W2 - O1O2). 


(3.1d) 



These particular quadratic expressions in Oj and o] are chosen because the commutation 



relations have a suggestive simple form: 

\Jri Ji 
[Q:Jr 

[Jr, Kg or L, 



0; 

ierst{Kt or L^), 

[Lr, Ls] = —i^rstJti 



(3.2a) 



(3.2b) 



(3.2c) 



Comparing and combining the above two sets of equations we see that the four operators 
Q and Jr are the photon-number conserving generators of U{2) (being respectively the 
generators of the U{1) and SU{2) parts of U{2))] while the six operators K^. and are the 
noncompact generators of squeezing transformations. 

The connection to SO {3, 2) becomes evident by regarding these generators as the various 
components of an antisymmetric set Mab — —Mba, where the indices A, B go over the range 
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(3.3) 



Then the commutation relations (|3.2D take on the de Sitter form 



[Mab, Mcd] = iigAcMsD - QbcMad + gAoMcB - QbdMca), 



(7AB = diag.(+l,+l, +1,-1,-1). 



(3.4) 



In this picture, Q, Jr are generators of rotations in the 4-5 plane and in the 1-2-3 subspace 
respectively, together making up the maximal compact 5*0(2) x 5*0(3) subgroup of 5*0(3, 2). 
With respect to the Fock or two mode photon number basis \ni,n2) for 7i the situation is 
that all these vectors with a fixed total number n = rii + n2, {n + 1) in all, are eigenvectors 
of Q with a common eigenvalue (n + l)/2; and they simultaneously provide the spin n/2 
representation of the 5*0(3) subgroup generated by Jr- On the other hand, the Kr and 
Lr,r = 1,2,3 are noncompact Lorentz boost generators in the r — 4 and r — 5 planes 
respectively. 

With this algebraic background, we go on to consider two-mode squeezing transforma- 
tions. It is a well known fact ||15| [T^ that each matrix 5 G 5p(4, 9ft) can be decomposed, 
globally and uniquely, into the product of two particular kinds of 5*p(4, 3ft) matrices: one 
factor belongs to the subgroup K, the other to a subset 11 defined in the following way: 



We shall hereafter denote elements in 11 by P, P', .... The decomposition mentioned above 
is then: 



n = {5 G 5p(4, gft) I 5^ = 5 = positive definite} C Sp{A, 3ft). 



(3.5) 



S G Sp{4, 3ft) 



S = P S{X,Y) 



p G n. 



X 



iY G U{2) 



(3.6) 
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with both factors being uniquely determined by S. This is the form in the present context 
of the famihar polar decomposition formula |]T^ for a general matrix. For the metaplectic 
operators U{S) we have the corresponding statement 

w(5)=w(P)w(5(x,r)), 

W(P) = exp i(real linear combination of K and L) , 
U{S{X,Y)) = exp i(real linear combination of Q and J) . (3-7) 

We may now identify precisely the most general squeezing transformation within the 
5*^(4, 9ft) framework, as the operator U{P) characterised by two numerical three-dimensional 
vectors k, I appearing as coefficients of K and L in the exponent: 

U{k, I) = exp [i{k -K + l- L)] . (3.8) 

Thus we reserve the name squeezing transformations for these elements of 11 within 5*^(4, 3ft), 
represented in the metaplectic representation by a single exponential factor. (It is well to 
keep in mind that 11 is not a subgroup, so the product of two such single exponential 
squeezing transformations is in general not a similar single exponential. This is analogous 
to the well known fact that the product of two 5*0(3, 1) Lorentz boosts is not just a boost 



but a boost followed or preceded by a rotation called the Wigner rotation |[18|). 

We may relate the decomposition ( |3.7|) of a general metaplectic transformation to a 
general quadratic Hamiltonian quite directly. Any such Hamiltonian containing both photon 
conserving and nonconserving terms with possibly time dependent coefficients would lead via 
the Schrodinger equation to a unitary finite time evolution operator which can be uniquely 
decomposed into the product form ( |3.7|) . Thus, integration of the Schrodinger equation leads 
in general to a specific passive factor and another specific squeezing transformation. In case 
the Hamiltonian is time independent and a combination only of the generators K and L, 
this evolution operator is already of the form U{P). 

Since we have a f/(2)-invariant squeezing criterion, as we have seen, elements of U{2) 
have no effect on the squeezed or nonsqueezed status of any given state. This means that 
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the U{2) transform, by conjugation, of a squeezing transformation is another squeezing 
transformation which should be regarded as equivalent to the first one. It is clear that, in 
any case, any equivalence relation among squeezing transformations as defined by us above 
should be based on processes which take one squeezing transformation to another. 

Now from the commutation relations ( ^.2b| ) we can see that the generators K, L, and the 
squeezing transformations U{k,l) defined in eq. (|3.8|) , behave as follows under conjugation 
by elements of U{2): 



e^" ■ '^{Kr or Lr) e~^" ' Rsr{a) {K^ or L 



cos 6 — sin 6 
sin 6 cos 



J 



k'l or I'l. = Rrs{c() {kg or I 



Rrsia) = Srs cos a + a^a 



1 — COS a 



+ ^rstOif 



sma 



a 



(3.9a) 



(3.9b) 



a = \a\ 



(3.9c) 



(3.9d) 



Here we have listed separately the effects of U{1) and SU{2) within U{2) on the generators 
and the squeezing transformations. The questions raised at the start of this Section can now 
be posed more precisely: If the set of squeezing transformations U{k, I) is separated into 
distinct nonoverlapping equivalence classes based on the U{2) action p.9| ), how can we con- 
veniently choose U (2)-invariant parameters to label these classes, and then pick a convenient 
representative element from each class? We answer these questions in this sequence. 

It is clear that we need to construct a complete set of independent expressions in k and /, 
invariant under both U{1) and SU{2) actions (|3.9b| , p.9c|) . We begin by defining the matrix 
M of scalar products among k and /, which is then SU{2) invariant: 
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M{k,l) 



^ k-k k-A 



V 



(3.10) 



k-l l-l 

This is a real, symmetric positive semi-definite matrix. With respect to U{1) action, we see 
from eq. ( |3.9b|) that M(fc, I) undergoes a similarity transformation by the rotation matrix 
of angle 6: 



M{k',f) = R{e) M{k,f) R{e)-\ 

COS 9 — sin 9 
sin 9 cos 9 

One now sees that there are two independent U (2) invariants that can be formed: 



R{9) 



( 

V 



(3.11) 



S5i(fc, I) = det M(fc, f)= \kA l\\ 
%{kj) = TiM{kJ)=\k\^ + \l\^; 



(3.12) 



and it is easily checked that there are no other invariants independent of these. 

Next let us tackle the problems of finding convenient parameters and representative 
squeezing transformations for the U (2) equivalence classes, one for each class. We see from 
eq. (|3.12|) that if > (i.e. 53i ^ 0) then the two vectors k and / are both nonzero 
and nonparallel; while if 53i = they are parallel (and one of them could vanish). These 
are therefore clearly different geometrical situations. Starting with the matrix M{k,l), we 
see from its U{1) transformation law ( |3.11|) that by a suitable choice of the angle 9 we 
can arrange the transformed matrix M{k', I') to be diagonal, and in the case of unequal 
eigenvalues to place the larger eigenvalue in the first position. This means that in each 
equivalence class of squeezing transformations there certainly are elements U{k, I) for which 
k ■ I = and \k\ > This still leaves us the freedom of action by SU{2). We may now 
exploit this freedom to put the (mutually perpendicular) vectors k and / into a convenient 
geometrical configuration. A look at the forms of the noncompact generators K and L in 
eq. (|3.1| ) suggests that we choose k and / as follows: 



fc = (0,a,0), / = (6,0,0), a>b. 



(3.13) 
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(A further reason for this choice will emerge shortly). 3i and can now be evaluated in 
terms of a, b to obtain the relations 

'^i{k,T) = a^b\ 
%ikj) = a^ + b^, 
a>b>0 , {a,b)^ (0,0). (3.14) 

With this parametrisation we can now say: there is a two-fold infinity of distinct equiva- 
lence classes of squeezing transformations for two-mode systems, each class corresponding 
uniquely and unambiguously to a point (a, b) in the octant a > 6 > in the real a — b plane, 
excluding the origin. Different points in the octant correspond to intrinsically distinct equiv- 
alence classes. Within an equivalence class determined by a point (a, 6), of course, one can 
connect different squeezing transformations U{k, l) by conjugation with suitable U{2) ele- 
ments. Given a squeezing transformation P G 11 C 5*^(4, 3?) its class (a, b) is determined by 
solving the equations 



/ „^ r . (a — b) (a + b). 

Tr (P) = 2 [cosh ^ ^ ^ + cosh ^ ^ \ 

Tr (p2) = 2[cosh(a - 6) + cosh(a + b)] 



(3.15) 



subject to the conditions on a and b appearing in eq. (|3.14|) . 

Then we have the following convenient two-mode squeezing transformation representing 
the equivalence class (a, 6): 

W(°)(«,^) =W^°Ha,0)W(°)(0,6), 
W(°)(a,0) =exp{iaK2) 



exp 



(giPi + ^2^2) 



W(o)(o,6) =exp(z6Li 
'ib 



exp 



(qiPi - q2P2) 



(3.16) 



The two factors U^'^\a, 0) and U^^\0, b) commute and may be written in either order, since 
according to eq. (|3.2c| ) the noncompact generators K2 and Li commute. 
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Finally one can easily calculate the symplectic matrix S^^\a, b) G 5*^(4, 3?), corresponding 
to the metaplectic operator U^^\a, b), by using eq. (^]8|). The result is: 

W(°)(a, b)-^ e W(°)(a, b) = S(°)(a, b) ^, 
5(°)(a, b) = diag. (e('^-'')/^ e^^+'^/', e-^^-^^l\ . (3.17) 



Now we can clarify that the particular choice ( p. 131 ) was dictated by the desire to have 



5''-°'' (a, 6) diagonal. This element of S'p(4, 3?) describes independent reciprocal scalings of 
the standard quadrature components of each mode. This amounts to showing geometrically 
that it is possible to diagonalise every P G 11 using conjugation by f/(2). 

We illustrate our classification scheme of two-mode squeezing transformations by giv- 
ing two examples. The extensively studied Caves-Schumaker [|| transformation uses the 
operator 

U^^'^){z) = exp a\ a\ - 03) . (3.18) 

By appearance, this attempts to involve or entangle the two modes maximally. In our 
notation this squeezing transformation corresponds to the generator combination 

z a\ — z* ai a2 = i{k ■ K + I ■ L)., 
k = -2{0,0,lmz), 

1= 2(0,0, Re 2). (3.19) 

Thus the invariant parameters a and b have values 

a = 2\zl 6 = 0. (3.20) 

The Caves-Schumaker squeezing transformations, and their U{2) conjugates, all taken to- 
gether, form a one-parameter family or one-dimensional line, in the a — b octant. In that 
sense they are a set of measure zero. 

Another interesting case is a squeezing transformation that refers essentially to a single 
mode but masquerades as a two-mode transformation: 



U{z] a, (3) = exp z{a*a\ + P*a2)'^ — z*{aai + (3a2y 

\a\^ + W = l. (3.21) 

After some simple algebra we find 

a, /?) = exp i{k ■ K + I ■ L) , 

k + il = 2z {-i{a'^ - p*^), + 2ia'p') . (3.22) 

The associated invariants and parameters are 

3l = 16|2|^ % = 8\z\'^, 

a = b = 2\z\. (3.23) 

These equivalence classes thus he along the line a = 6 in the octant, again a one-parameter 
family of zero measure. Our results are depicted in Figure |l|. 

We note here that the size of an equivalence class depends sensitively upon the point 
(a, b). Since for a 7^ 0, 6 7^ and a b none of the generators of U (2) or a linear combination 
of them commute with aK2 + bLi, we have a full four parameter equivalence class. For the 
cases a 7^ 0,6 = 0, and a = b, respectively, the vanishing of the commutators [J2,K2] and 
[Q + J3, K2 + Li] leads to reduction of the dimensionality of the equivalence class to three. 

We conclude this Section with a few comments. Each point (a, b) in the octant denotes 
an equivalence class of squeezing transformations, whose dependences on a and b would be 
of physical significance and would show up in a variety of [/(2)-invariant properties. The 
two-mode transformations so far discussed in the literature lie basically along the two lines 
shown in Figure |l|. In this sense, most of the intrinsically distinct two-mode transformations, 
their effects on various states, etc., remain to be explored. Those equivalence classes (a, 6) 
for which a > b involve the two modes in an essential way. We may say purely qualitatively 
that the distance of the point (a, b) from the line a = 6, or perhaps better the expression 
(1 — b/a), is a measure of the extent to which two independent modes are involved in the 
transformation. In this sense, as remarked earlier, the Caves-Schumaker transformations 
involve two modes maximally. 

19 



IV. SQUEEZED COHERENT AND THERMAL STATES FOR TWO MODES 



The general two mode coherent state with complex two-component displacement a = 
(ai, 02) is defined by 

la) = exp (ji ■ a'' — a* ■ aj |0, 0) 

= exp - ■^|a2p^ exp (^aia| + a24) |0,0). (4.1) 

For this state the means of the quadrature components C,a do not vanish in general: 

= v^(Reai,Rea2,Imai,Ima2)^- (4.2) 

The variance matrix is however independent of a: 

V{\a)) = v{\0)) = ^U.,. (4.3) 

The most general squeezed coherent state is obtained by applying U{P) for some P G 
n C 5*^(4, 3?) to \a) for some a. This ly({P) is conjugate, via some f/(2) element, toW(°)(a,6) 
for some a, b. Now the effect of a U{2) transformation on \a) is to give us another coherent 
state \<y'), dl being the f/(2) transform of a. But the variance matrix is in any case a- 
independent. To examine the [/ (2) -invariant squeezing condition, therefore, it suffices to 
examine the particular class of squeezed coherent states 

|a;a,6) = W(°)(a,6) (4.4) 



From eqs ( p. 12 , 3. 17 ), the calculation of the variance matrix for this state is trivial, and it is 



in fact diagonal: 

y(|a;a,6)) = S^''\a,h)V {\a)) S^^\a,b) 
= ^S^'\2a,2b) 

= i diag. (e^'^-^), e('^+^), e^''-'^), e-("+^)) . (4.5) 

Since a and b are nonnegative, and in addition a -|- 6 > 0, we see that the least eigenvalue of 
this variance matrix is 
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(V(|a;a,6))) = le-('^+^)<i 



(4.6) 



These states are thus always squeezed. 

If we apply any passive U{2) transformation S{X,Y) to any one of the states \a;a,b) 
defined above, the variance matrix will in general change as V ^ V = S{X, Y) V S{X, Y)'^; 
but its eigenvalue spectrum, and in particular i{V), remains unaltered. Thus all the states 
symbolically written as U{U (2)) \a; a, b), for various U (2) elements, are squeezed to the same 
extent as \a; a, b), and have i{V) given by eq. ( |4.6|) . 



The Schrodinger wave functions for the subfamily of squeezed coherent states ( [4.4|) are 
particularly simple, since they are products of single mode squeezed coherent state wave- 
functions: 



q2\a; a, b) = ip''°\q[; ai, a -b) ip^^\q2, a2, a + b) 



(0)l 



^W(g';a,a) 



-a/4 



1/4 



exp 



ia Ima — - {q e "'^^ — \/2a 



(4.7) 



(For a general state U{U{2)) |a;a,6), we do not expect such a product form). When we 
set 6 = (Caves-Schumaker limit), both factors show the same amount of squeezing; while 
when we set a = 6 (essentially single mode situation) we see squeezing only in the factor 
referring to the second mode. These features are as we would have expected. 

The next example we look at is the case of a two-mode thermal state subjected to 
squeezing. The motivation in making this choice is that the starting density operator is 
explicitly U{2) invariant. The normalized density operator corresponding to inverse tem- 
perature (3 = huj/kT is described in the Fock and Schrodinger representations by: 



Po(/5) 



Po(gi,g2)9i,?2;/^) 



[I — e exp [— /3(ai ai + 02 ^'2] 
00 



(4.8a) 



tanh^ 



-\ (^tanh ^ + coth ^ {q'Y + + gf + gf ) 
(^tanh^-coth|^ {qW; + qWD 



(4.8b) 
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with U (2) invariance expressed by 

e'^Q Po{f3) e-'^Q = e^" " ^ po{f3) e''^ " J = po{f3). (4.9) 

Therefore it suffices to examine the properties of the density operator obtained by conju- 
gating poi(3) with W(°)(a, b): 

p{p;a,b) =U'^^\a,b) po{f3) U^^\a,b)-\ (4.10) 

The most general squeezed thermal state is evidently 

W(f/(2))p(/5;a,6)W(f/(2))-\ (4.11) 

but this has the same squeezing properties as p(/5; a, b). 



For the thermal state Pq{(3) the variance matrix is well known |p!0| : 



V{pm) = \ COth^l4x4 (4.12) 

Therefore for the particular set of squeezed thermal states ( [l.l(j| ), we have diagonal variance 
matrices: 



V{p{(3-a,b)) =5(°)(a,6)\/(po(/5)) S^'\a,b) 
= ^ coth^ 5(°)(2a,26) 



2 2 
The least eigenvalue is evidently 



]- coth I diag. (e^"^-') , e('^+') , e^'"'^) , e"'^"') . (4. 13) 



£(\/) = 1 coth^ e-('^+'\ (4.14) 

so for a given temperature, squeezing sets in when 

a + 6>lncoth^. (4.15) 

In Figure ^ this region consists of all points in the a — b octant to the right of the line 
a + b = lncoth/?/2, which is a line perpendicular to the line a = b and at a distance 
In coth(/?/2) /V2 from the origm. 
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V. DETECTION SCHEMES AND ROLE OF U{2) TRANSFORMATIONS 



We have so far not specified in any detail the two orthogonal modes of radiation being 
subjected to squeezing. Let us at this point consider a situation well studied experimentally 



by the heterodyne detection scheme [jl9l . Here the two modes differ only slightly in frequency, 
but are otherwise similar. In this kind of experimental arrangement, what is actually mea- 
sured is the fluctuation of a certain photocurrent, and this in turn gives the fluctuation or 
variance of the g-quadrature component of a particular (passive) combination of the original 
modes. The combinations of ai and 02 that are involved form the one-parameter family 



a{^) = ^ («i + «2) e-*^/2 



< < 47r. 



(5.1) 



This can be regarded as the first component of a [/(2)-transformed pair a[, a'g: 





1 ( 







Q-ii>/2 g-jV/2 



ai 



a{ip). 



(5.2) 



The hermitian quadrature component whose fluctuation is measured is 



(5.3) 



The only experimentally adjustable parameter here is the angle ip. The family of U{2) 
elements realised in the heterodyne scheme is thus only the one-parameter set given in eq. 
( ^.2|) parametrised by ip, and belonging to SU{2): 



1 

V2 



( g-#/2 g-iV/2 \ 



_g#/2 gjV/2 



G SU{2). 



(5.4) 



We notice that this is not a one-parameter subgroup of S\J (2); in particular even the identity 
element of the group is not contained here. 

With this description of the heterodyne setup in our framework, let us see to what extent 
it can be used to detect t/ (2) -invariant squeezing. Now a general two-mode state p with 
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variance matrix V, even if it is squeezed in the intrinsic sense of eq. ( p.l7|) , may not be 



manifestly squeezed. That is, it may happen that Vaa > 1/2 for all a = 1, . . . , 4. As our 
discussion in Section (|I|) shows, we need to be able to experimentally realise a general U (2)- 
transformation applied to the state p, and change its variance matrix to a form where one 
of its diagonal entries (say the leading one) becomes less than 1/2. However the heterodyne 
method is generally unable to do this job for us, as it can only realise the one-parameter 
subset of SU{2) transformations Uni'ip) for < ip < An. 

In the two examples of squeezed coherent states and squeezed thermal states studied in 
the previous Section, we have a family of states related to each other by conjugation with 
U{2) for each point in the a — b plane. Each equivalence class has appropriate dimensionality 
depending upon the point (a, b) as explained in Section |T|. It turns out that for each (a, b) 
the heterodyning scheme can detect squeezing in only a one parameter subset of the family 
of states. Although heterodyning detection covers the whole a — b plane it does not reach all 
the states corresponding to each point in the a — b plane. For example, in the representative 
chosen in eqs ( |4.4| , ^.1CI| ), for which the variance matrix is already diagonal, the squeezing 
cannot be detected by this scheme because of the absence of the identity in Uh^. It should 
be possible to detect squeezing in these states by a suitably modified scheme. We wish to 
emphasize that there is a definite need to be able to experimentally implement the most 
general elementof U{2). This would allow the experimenter to detect the degree of squeezing 
unambiguously, if the state is squeezed, without any prior knowledge of the elements of the 
initial variance matrix. 

Having stressed the need to implement arbitrary U{2) transformations on the two modes 
of radiation in order to reach the proper quadrature to exhibit squeezing, we now describe 
how it can be achieved in some situations. We discuss two particular cases of the two modes 
involved, the first when the two modes have the same frequency but different directions 
of propagation, and the second when the modes have the same frequency and direction of 
propagation but different polarisations. 

The experimental setup for the first case is shown in Figure We achieve an arbitrary 
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U{2) transformation on the two modes by using a Mach-Zehnder interferometer with two 



50:50 beam sphtters (BSi and BS2) and appropriate phase shifters The input modes 
with annihilation operators ai and 02 are subjected to equal and opposite phase shifts by 
angles and —(f), then the modes are mixed in the beam splitter BSi, the mixed modes again 
undergo equal and opposite phase shifts by angles 6 and —6 and further mixing through 
the beam splitter BS2. Finally they undergo unequal phase shifts by angles ipi and 4'2- If 
the annihilation operators at the output are a[ and a'g then all the above operations when 
combined are implemented through the transformation 



\(^2 J 



eiini^i) cos e -ie-*(<^-'^i) sin 9 



ya2y 



(5.5) 



relating the two sets of annihilation operators. The above matrix is the most general U{2) 
transformation matrix. We note here that if 1^2 = —^1 then the transformation matrix is the 
most general SU{2) transformation. Going from SU{2) to U{2) is just a matter of overall 
phase and can also be achieved by free propagation. 

For the second case when the two modes differ only in polarisation we achieve the arbi- 
trary U{2) transformation by using two quarter wave plates (Qi, Q2) and a half wave plate 
H as shown in Figure |^. The detailed discussion of this set up is given in ||21]]. It turns out 
that the configuration Q-H-Q is not the only one but Q-Q-H and H-Q-Q also accomplish the 
same result, as shown in [^. We basically have three elements, a quarter wave plate Qi, a 
half wave plate H and a quarter wave plate Q2 all three of them being coaxially mounted 
and with their slow axes in the x — y plane making angles of a, (3 and 7 respectively with 
the X axis. The two modes having annihilation operators ai and 02 moving along the z 
direction pass through this arrangement. If the annihilation operators at the output are 
a'l and a'2 then they are related to the operators at the input by an SU{2) transformation 
given in terms of a, (3 and 7. By changing these parameters one can reproduce any desired 
SU{2) element. As has been pointed out earlier going to U{2) now is just a matter of free 
propagation. 
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In both the above cases, by going to the proper U{2) element we can make the squeezing(if 
present) manifest and bring it to the leading diagonal element of the variance matrix i.e. in 
the quadrature q[ = + The squeezing in this quadrature can now be measured 

by any standard one mode detection method. 

These remarks show on the one hand the way the heterodyning scheme fits into our 
general analysis, and on the other hand the need to devise new schemes capable of realising 
all elements of U{2), tailored to the definition of the two modes involved. 

VI. CONCLUDING REMARKS 

We have presented a classification scheme for two-mode squeezing transformations, based 
on the structure of the real four dimensional symplectic group 5*^(4, 9ft), and the separation 
of its elements into passive (compact) and active (noncompact) types. The structure and 
action of the maximal compact subgroup U{2) in 5*^(4, 9ft), and the ?7(r;,)-invariant squeezing 
criterion formulated elsewhere for a general n-mode system, have guided our considerations. 
All our work is in the metaplectic unitary representation of S'p(4, 3ft); and the local isomor- 
phism 5*^(4, 9ft) ~ 50(3, 2) has led to a convenient geometric picturization of the situation. 

As emphasized in Section |T| the squeezing transformations U{P), P G 11 C 5*^(4, 9ft), 
do not form a subgroup of 5'p(4, 3ft). The breakup of these transformations into equivalence 
classes, based on the effect of conjugation by elements of U (2) is the only natural available 
classification procedure. This is because the definition of equivalence classes for any set of 
objects has to be based on an equivalence relation defined on that set. Thus, we have treated 
two elements P, P' G 11 as intrinsically equivalent if 

P' = S{X,Y) P S{X,Yf for some X-iYeU{2). (6.1) 

It should however be realised that the detailed effects of action by W(P) and U{P') on a 
general initial two mode state po? as seen in the changes caused in the variance matrix V{pq), 
need not be identical. Since this is a subtle and important point we spell it out in detail. 
Starting from a general state po? action by a squeezing transformation leads to a new state 
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p = U{P) p^U{P)-\ (6.2) 
As seen in Section |T|, any hi{P) is expressible in terms of a representative element U^^\a, h) 



as 



U{P) = U{S{X, Y)) U^^\a, b) U{S{X, Y))-^ (6.3) 
for suitable X — iY G U{2). Therefore we have 

p = U{S{X, Y)) U^^\a, h) U{S{X, F))"^ po U{S{X, Y)) U^^\a, b)-^ U{S{X, Y))-\ (6.4) 
This leads by eq. ( |2.12|) to the relation 



V{p) = S{X, Y) 5(°)(a, b) S{X, Y f V{po) S{X, Y) S^^\a, b) S{X, Y f (6.5) 

between the two variance matrices. Now the first and last factors on the right hand side 
here have no influence on the spectrum, and so on the least eigenvalue, of V{p). Therefore 
the squeezed or nonsqueezed nature of p is actually determined by the least eigenvalue of 
the matrix 

5(X, Yf V{p) S{X, Y) = 5(°)(a, b) 5(X, Yf V{po) S{X, Y) S^'^\a, b). (6.6) 

But now the right hand side is in general dependent not only on the invariant parameters 
a, b but also on X, Y . In the examples studied in Section namely where po is a coherent 
state or an isotropic thermal state, V{pq) happens to be a multiple of the identity matrix, 
so that on the right hand side of eq. (|6.6|) the dependence on X, Y cancels. But this need 
not happen in general. Thus for instance if we take for po an anisotropic thermal state 
with unequal temperatures for the two modes, we have only U{1) x [/(I), rather than f/(2), 
invariance for this po; so the least eigenvalue ^{y{p)) of V{p) will depend on a, b and on two 
out of the four U{2) parameters present in X — iY . One can easily convince oneself that 
the only situation where ^(X, Yf V{pq) S{X, Y) = V{po) independent of X and Y is when 
V{po) is a multiple of the unit matrix; and the isotropic thermal states do reproduce all such 
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cases. Therefore a more detailed study of the effect of a general squeezing transformation 
on initial states po with nontrivial V{pq) is of considerable interest. 

A related important question is the following: Let us take two squeezing transformations 
Pi, P2 G n belonging to equivalence classes (ai,6i), (02,^2) respectively, which could coin- 
cide. The product P1P2 will in general be of the form S{X, Y) P with S{X, Y) & U (2) and 
P G n. If P belongs to the equivalence class (a, b) we wish to determine this class in terms 
of Pi and P2. Using the fact Tr{P^) = Tr((PiP2)(PiP2)^) and eq. (|3l5D we arrive at the 
following relations 

2[cosh {a-b) + cosh (a + b)] = Tr{{PiP2){PiP2f) 
2[cosh2(a -b) + cosh2(a + b)] = Tr{{{PiP2){PiP2ff) (6.7) 

which can be solved to find (a, b). We note here that (a, b) do depend not only upon (ai, 61) 
and (02, 62) but also on the actual elements chosen from each of these classes. So we do not 
have a notion of class multiplication among these equivalence classes. 

Finally we call attention to our considerations in Section |V| and to the need for being 
able to experimentally implement or realize general passive elements of the subgroup U (2) 
of 5*^(4, 5R) for each given choice of the independent modes in a two-mode system. Once 
this is achieved, for any given state we can bring out in an explicit or manifest fashion 
its squeezing nature (provided it is squeezed) by altering its variance matrix and making 
the least eigenvalue appear in the leading position on the diagonal. This also means that 
we would be able to experimentally measure the fluctuation in the quadrature variable 
isolating the least eigenvalue. As extensively discussed elsewhere, these considerations which 
exploit the richness of the geometry underlying the symplectic group do not require complete 
diagonalization of the variance matrix at all |[T0[| . 

In the following paper we shall examine ?7 (2) -invariant properties of two mode squeezed 
states which go beyond the level of second order moments of the quadrature operators. 
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FIGURES 




FIG. 1. Equivalence classes of two-mode squeezing transformations, Caves-Schumaker (C-S) 
and single mode limits, squeezed thermal region. 
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FIG. 2. Mach-Zehnder Interferometer implementing arbitrary U{2) transformation on two 
modes at the same frequency but differing in their direction of propagation. BSi and BS2 are 
50:50 beam spHtters and thick hnes are phase shifters by angles indicated, ai, 02 are the an- 
nihilation operators at the input port and a[, are the annihilation operators at the output 
port 



34 



X 

y 

^ * z Qi H Q2 



a 









7 




— 


0102 











FIG. 3. Implementation of arbitrary SU{2) element on two modes with the same frequency 
and directions of propagation, but different polarisations. Here Qi, Q2 are quarter wave plates and 
H is the half wave plate, a, /3, and7 are the angles which the slow axes of Qi, H and Q2 make with 
the X axis respectively, ai, 02 are the annihilation operators at the input port and a'l, a'2 are the 
annihilation operators at the output port 



35 



